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"What  then,"  asked  St.  Augustine,  "is  time?  If  no  one  asks  me,  I  know 
what  it  is.  If  I  wish  to  explain  it  to  him  who  asks  me,  I  do  not  know." 

We  have  learned  a  few  things  since  St.  Augustine,  Einstein  has  taught  us  a 
lot.  But  still  there  are  a  lot  of  unanswered  questions.  In  particular, 
how  do  you  measure  time?  It  intrigues  me  that  we  never  measure  time;  we 
measure  time  differences,  i.e.  the  time  difference  between  two  clocks.  I 
know  of  no  way  to  measure  the  time  of>a  clock.  I  can  measure  the  time  of 
an  event  with  reference  to  a  particular  clock.  Another  intriguing 
question  is,  if  time  cannot  be  measured,  is  it  physical  or  is  it  an 
artifact?  We  conceptualize  some  of  the  laws  of  physics  with  time  as  the 
independent  variable.  We  attempt  to  approximate  our  conceptualized  ideal 
time  by  inverting  these  laws  so  that  time  is  the  dependent  variable.  The 
fact  is  that  time,  as  we  now  generate  it,  is  dependent  upon  defined 
origins,  a  defined  resonance  in  the  cesium  atom,  interrogating  electron¬ 
ics,  induced  biases,  and  random  perturbations  from  the  ideal.  Hence,  at  a 
significant  level,  time  —  as  man  generates  it  by  the  best  means  available 
to  him  —  is  an  artifact.  Corollaries  to  this  are  that  every  clock 
disagrees  with  every  other  clock  essentially  always,  and  no  clock  keeps 
ideal  or  "true"  time  except  as  we  may  choose  to  define  it.  Frequency  or 
time  interval,  on  the  other  hand,  is  fundamental  to  nature;  hence,  the 
definition  of  the  second  can  approach  the  ideal.  Noise  in  nature  is  also 
fundamental.  Characterizing  the  random  variations  of  a  clock  opens  the 
door  to  optimum  estimation  of  environmental  influences  and  to  the  design 
of  optimum  combining  algorithms. 


45 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

DEC  1985 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-1985  to  00-00-1985 


4.  TITLE  AND  SUBTITLE 

Characterization,  Optimum  Estimation,  and  Time  Prediction  of  Precision 
Clocks 

6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

National  Bureau  of  Standards, Boulder, CO, 80303 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 


13.  SUPPLEMENTARY  NOTES 

Proceedings  of  the  Seventeenth  Annual  Precise  Time  and  Time  Interval  (PTTI)  Applications  and  Planning 
Meeting,  Washington,  DC,  3-5  Dec  1985 


14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


17.  LIMITATION  OF 

18.  NUMBER 

ABSTRACT 

OF  PAGES 

Same  as 

24 

Report  (SAR) 

19a.  NAME  OF 
RESPONSIBLE  PERSON 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Let's  define  some  terms^-^:  Beginning  with  a  sine  wave  voltage  with 
frequency  v(t)  as  the  variable  frequency  output  of  a  precision  oscillator, 
we  may  write  V(t)  =  VQ  sin(2ir  v(t)*t),  where  we  assume  that  amplitude 
fluctuations  are  negligible  around  VQ,  and  where  v(t)  is  the  average 
frequency  from  0  to  t.  We  can  redefine  this  equation  with  being  a 
constant  nominal  frequency  and  place  all  of  the  deviations  in  the  phase, 
V(t)  =  VQ  sin  ( t  +  $(t)).  We  then  define  a  quantity 
y(t)  -  (v(t)  -  which  is  dimensionless  and  which  is  the  fractional 

frequency  deviation  of  v(t)  from  the  nominal  value.  We  can  integrate 
y(t)  to  get  the  time  deviation,  x(t),  which  has  the  dimensions  of  time. 
From  the  above  we  can  write  the  time  deviation  of  a  clock  as  a  function  of 
the  phase  deviation:  x(t)  ■  4>( t) /2irv0. 


Why  do  we  have  time  deviations?  We  conceptualize  two  categories: 
systematica,  such  as  frequency  drift  (D),  frequency  offset  (yQ)  and  time 
offset  (xQ);  and  then  random  deviations  e(t),  which  are  not  considered  to 
be  deterministic: 


x(t)  -  xG  +  yQ»t  +  1/2  Dt2  +  e(t).  (1) 

Note,  the  quadratic  in  the  D  term  is  because  x(t)  is  the  integral  of  y(t), 
the  fractional  frequency.  In  Figure  1  we  have  simulated  two  systematic 
cases:  one  a  clock  with  frequency  offset,  and  another  case  with  a 
negative  frequency  drift,  systematlcs  caused  by  environmental  influences 
are  also  very  important.  Figures  2  through  6  summarize  some  of  the 
important  systematic  influences  on  precision  oscillators.  An  important 
set  of  systematic  deviations  are  modulation  side  bands,  e.g.  60  Hz, 

120  Hz,  daily  and  annual  dependences,  those  induced  by  vibrations,  etc. 


The  random  deviations  of  precision  oscillators  can  typically  be  catego¬ 
rized  by  power  law  spectra,  Sy(f)  -  fa  where  f  is  the  Fourier  frequency 
and  a  takes  on  integer  values,  i.e.  -2,  -1,  0,  1,  2^’2,3»2^.  Figure  7 
shows  noise  samples  corresponding  to  these  different  power  law  spectra  and 
Table  1  shows  the  nominal  range  of  applicability  of  these  power  law 
models.  Given  a  time  deviation  plot  x(t)  for  the  time  difference  between 
a  pair  of  clocks  or  a  clock  against  some  primary  reference,  and  some 
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sample  time  t  (see  Figure  8  for  example),  the  average  fractional  frequency 
for  each  interval  is  the  time  difference  at  the  end  of  the  interval  minus 
that  at  the  beginning  of  the  interval  divided  by  t.  We  can  thus  construct 
a  set  of  discrete  frequency  values  over  such  a  data  set  from  this  time 
deviation  plot.  We  can  calculate  a  classical  standard  deviation  for  these 
values,  but  one  can  show  that  for  some  kinds  of  power  law  spectra 
encountered  in  precision  oscillators  that  the  standard  deviation  is 
divergent-"’-! ,  i.e.,  it  does  not  converge  to  a  well  defined  value,  and  it 
is  a  function  of  data  length^-2-!.  Hence,  the  standard  deviation  should  not 
be  used  for  characterizing  clocks.  An  IEEE  subcommittee  has  recommended 
Sy(f)  in  the  frequency  domain  and  a  measure,  Oy2(i)  in  the  time  domain^-1-!. 
The  latter  has  come  to  be  called  the  two-sample  variance  or  the  Allan 
variance.  The  convergence  of  Oy(T)  has  been  ver if ied^ for 
the  power  law  spectra  of  interest  in  precision  oscillators.  It  is  defined 
as  follows^-1  -! : 

oy2(T)  -  j  <Ay)2>  (2) 

where  Ay  is  the  difference  between  adjacent  fractional  frequency 
measurements  each  sampled  over  an  interval  t  and  the  brackets  <> 
indicate  the  infinite  time  average  or  expectation  value. 


A  pictorial  description  is  shown  in  Figure  9  for  a  finite  data  set.  A 

data  set  of  the  order  of  100  points  is  quite  adequate  for  convergence  of 

0y(O,  though  of  course  the  confidence  of  the  estimate  will  typically 
improve  as  the  data  length  increases^-!.  Given  a  discrete  set  of  stored 
data,  the  value  of  t  can  be  varied  in  the  software^-^ .  If  tq  is  the  data 
spacing  for  the  stored  data  set,  y^,  from  the  measurement  system,  one  can 

change  to  t  =  nx0  by  averaging  n  adjacent  values  of  yjto  obtain  a  new 

fractional  frequency  estimate  with  sample  time  t  as  input  to  Equation  (2). 
Hence,  in  a  very  convenient  way  one  can  calculate  0y(x)  as  a  function  of 
x,  which  will  be  shown  to  be  very  useful  in  a  moment.  For  a  finite  data 
set,  Equation  (2)  then  becomes 
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0y2(O 


2t  (M+1-2n) 


M+1-2n 

l 


+  x,V 


where  the  yk  are  average  fractional  frequencies  averaged  over  t  «  m0, 
starting  at  k,  and  the  xks  are  the  discrete  time  deviation  measurements, 
k  _ 

xk  -  kx0  l  yr 

i  =  1 

Equation  (3)  is  obtained  from  a  first  difference  on  frequency,  and 
equation  (4),  from  the  second  difference  on  the  time;  they  are  mathemati¬ 
cally  identical,  which  gives  us  the  option  of  using  frequency  or  time 
data. 

For  power  law  spectra,  y  is  constant  for  a  particular  value  of  a,  where 
Oy2(-r)  -  The  relationship  of  the  spectral  density,  Sy(f)  -  fa, 

exponent  a  versus  p  is  p  =  -a-1  (-2<  a  £1)  and  y  =  -2  (ail). 

For  example,  oy(T)  is  proportional  to  t  (which  is  typical)  for  cesium, 
rubidium  and  passive  hydrogen;  then  y  has  the  value  of  -1,  and  hence  a  has 
the  value  of  0  (white  noise  frequency  modulation) .  This  is  the  classical 
noise  exhibited  by  most  atomic  clocks  for  t's  beyond  a  few  seconds  and  in 
this  case  Oy(x0)  is  equal  to  the  classical  standard  deviation.  Fortu¬ 
nately  for  most  cases  where  Til  second  the  relationship  y  =  -a-1  is 
applicable. 

We  have  an  ambiguity  at  y  =  -2;  we  cannot  tell  whether  we  have  flicker 
noise  phase  modulation  (PM)  or  white  noise  PM.  We  can  avoid  this  problem 
by  realizing  that  in  this  region  oy(t)  depends  on  the  measurement 
bandwidth^-2  .  One  can  construct  a  variable  software  bandwidth,  f  ,  by 
realizing  the  following^-®  .  In  any  measurement  system  a  hardware 
bandwidth,  f^,  exits  through  which  we  measure  the  phase  or  the  time 
difference  between  a  pair  of  clocks;  define  =  1/2-irf^.  In  other  words 
is  the  sample  time  period  through  which  we  sample  the  data.  If  we 
average  n  time  or  phase  readings,  we  increase  the  time  period  to  =  ts. 
But  ts  *■  l/2irf s,  where  fs  =  f^/n,  i.e.,  we  narrow  the  effective  bandwidth, 
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f  ,  in  the  software  by  n.  In  other  words  fg  =  f^/n  gets  smaller  as  we 
average  more  values;  i.e.  increase  n  (t  =  m0).  We  can  therefore 
construct  a  second  difference  composed  of  time  deviations  so  averaged,  dnd 
then  define  a  modified  0y2(x)  that  will  remove  the  ambiguity  through 
bandwidth  variation: 


9  1  N-3n+l  /n+j 

Mod.  o  ( x )  * — -  I  (  £ 

y  2x  n  (N-3n+l)  j-1  \  i-J 


(xi+2n  "2xi+n  * 


(5) 


where  N  *  M+l,  the  number  of  time  deviation  measurements  available  from 
the  data  set.  And  now  if  Mod.0y2(t)  -  iy,  then  y'«  -a-1  ( 1£  a  <3)^’10-^ 


We  typically  employ  Mod.Oy(t)  as  a  subroutine  to  remove  the  ambiguity  if 
Oy(i)  ~  t” ^  because  the  y'  dependence  only  approximates  that  given  by 
equations  (3  and  4)  for  o  <1.  But  for  a  =  2  and  1,  y  /2  exactly  equals 
-3/2  and  -1  respectively,  providing  a  clean  differentiation  between  white 
noise  PM  and  flicker  noise  PM. 


Table  2  illustrates  why  one  should  not  use  the  classical  standard 
deviation  to  characterize  clocks.  For  the  different  kinds  of  noise 
processes  we  list  the  classical  standard  deviation  of  the  time  deviations 
and  the  classical  standard  deviation  of  the  fractional  frequency 
deviations  as  a  function  of  Oy(-r)  (the  square  root  of  the  Allan  variance). 
The  divergent  nature  of  either  classical  standard  deviation  is  apparent, 
and  even  for  classical  white  noise  FM  the  standard  deviation  is  apparent, 
and  even  for  classical  white  noise  FM  the  standard  deviation  of  the  time 
diverges  as  the  square  root  of  the  data  length  i.e.  the  number  of  samples 
N[2]. 

Using  Oy(t)  or  Mod.Oy(i)  we  can  characterize  typical  power  law  processes. 
We  then  have  the  opportunity  of  determining  optimum  estimates  of 
values  by  employing  the  statistical  theorem  that  the  optimum  estimate  of  a 
white  noise  process  is  just  the  simple  mean. 
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For  example,  consider  the  very  common  and  very  important  case  of  white 
noise  FM  typically  found  on  the  signals  from  cesium  standards,  rubidium 
standards  and  passive  hydrogen  masers.  The  optimum  estimate  of  the 
frequency  is  the  simple  mean  frequency,  which  is  equivalent  to 
(xN  -x^/MTq.  It  is  still  all  too  common  within  our  discipline  to  see  our 
colleagues  erroneously  determining  the  frequency  for  these  kinds  of  , 

oscillators  by  calculating  the  slope  from  a  linear  leasts  squares  fit  to 
the  time  deviations  and  quoting  the  standard  deviation  around  that  fit  as 
a  measure  of  the  clock  performance.  There  are  three  problems  in 
proceeding  this  way.  First,  the  frequency  estimate  is  not  optimum  In  a  j 
mean  square  error  sense  and  is  equivalent  to  throwing  away  about  20$  of 
the  data,  Increasing  the  cost  in  the  case  of  a  calibration.  Second, 
the  standard  deviation  diverges  as  the  square  root  of  the  data  length  and 
third,  the  standard  is  significantly  dependent  on  the  filter  (linear 
leasts  squares)  as  well  as  the  clock  deviations.  On  the  other  hand  such  a 
filter  can  be  useful  for  assessing  outliers.  The  optimum  "end  point" 
method  outlined  above  has  the  risk  that  if  either  of  the  points  is 
abnormal,  i.e.  the  model  fails,  the  result  will,  of  course,  be  adversely 
effected,  so  such  a  filter  is  useful  to  assess  whether  there  are,  outliers 
—  paying  especial  attention  to  the  end  points. 

There  are  other  useful,  and  maybe  not  so  obvious,  optimum  estimators  at 
the  conclusion  of  a  data  set:  (1)  Given  white  noise  PM,  the  best  time 
error  estimate  is  the  simple  mean  of  the  time  deviations,  the  frequency 
estimate,  then,  is  the  slope  from  a  linear  least  squares  fit  to  the  time 
deviations,  and  the  frequency  drift,  D,  is  determined  from  a  quadratic 
least  squares  fit  to  the  time  deviations  per  equation  (1).  (2)  Given 

white  noise  FM,  the  optimum  estimate  of  the  time  is  the  last  value,  the 
optimum  frequency  estimate  is  outlined  in  the  previous  paragraph  and  the 
optimum  frequency  drift  estimate  is  from  a  linear  least  squares  fit  to  the 
frequency.  (3)  Given  random  walk  FM,  the  optimum  time  estimate  is  the 
last  value,  and  optimum  frequency  estimate  is  obtained  from  the  last  slope 
of  the  time  deviations,  and  the  optimum  frequency  estimate  is  from  the 
last  slope  from  the  time  deviations,  and  the  optimum  frequency  drift 
estimate  is  calculated  from  the  mean  second  difference  of  the  time 
deviations.  Caution  needs  to  be  exercised  here  for  typically  there  Will 
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be  higher  frequency  conponent  noise  in  a  real  data  stream,  such  as  white 
noise  FM,  along  with  the  random  walk,  and  these  can  significantly 
contaminate  the  drift  estimate  from  a  mean  second  difference.  If  random 
walk  FM  is  the  predominant  long-term,  power- law  process,  which  is  often 
the  case,  then  the  noise  can  be  handled  by  calculating  the  second 
difference  from  the  first,  middle  and  end  time  deviation  points  of  the 
data.  The  flicker  noise  cases  are  significantly  more  complicated,  though 
filters  can  be  designed  to  approximate  optimum  estimation  ^11  ,12,133^  jn 
the  limit  as  the  data  length  increases  without  limit,  the  time  is  not 
defined  for  flicker  noise  PM,  an  the  frequency  is  not  defined  for  flicker 
noise  FM.  This  has  some  philosophical  implications  for  the  definitions  of 
time  and  frequency  unless  some  low  frequency  cutoff  limits  exist.  If 
significant  frequency  drift  exists  in  the  data,  it  should  be  optimally 
subtracted  from  the  data  or  it  will  bias  the  long  term  values  of  Oy(x): 

°y(T)  *  •  (6) 

Once  the  power  law  spectra  are  deduced  for  a  pair  of  oscillators,  then  one 
can  also  develop  an  optimum  predictor.  Table  3  gives  both  the  optimum 
prediction  values  for  the.  various  relevant  pure  power  law  spectra,  as  well 
as,  their  assymtotic  forms.  Special  forecasting  techniques  must  be  used 
for  optimal  prediction  when  combinations  of  these  processes  are  present. 

To  illustrate  how  these  concepts  relate  to  real  devices,  Figure  10  shows 
a  °y(x)  diagram  for  some  interesting  state-of-the-art  oscillators,  and 
Figure  11  shows  the  rms  time  prediction  errors  for  the  same  set  of 
oscillators. 

In  conclusion  it  is  clear  that  classical  statistics  does  not  allow 
characterization  of  common  kinds  of  random  signal  variations  found  in 
precision  oscillators.  The  two-sample  or  Allan  variance  provides  an 
efficient  and  convergent  measure  of  the  power  law  spectral  density  models 
useful  in  characterizing  most  of  these  oscillators.  Once  characterized  we 
can  calculate  optimum  time  and  frequency  estimates  as  well  as  predicted 
values.  Characterizing  the  random  variations  also  provides  near  optimum 
estimation  of  systematic  effects,  which  often  cause  the  predominant  time 
and  frequency  deviations.  For  example,  if  we  wanted  to  optimally 
determine  the  temperature  dependence  with  the  temperature  set  at  two 
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different  values,  we  would  set  at  one  temperature  and  measure  the 
frequency  against  a  reference  for  a  time  im,  corresponding  to  the  t  for 
the  minimum  Oy(i)  value.  We  would  then  change  the  temperature  to  the 
other  value  and  repeat  the  measurement  with  the  same  criteria  and  note  the 
resultant  Ay  between  the  two  optimally  determined  frequency  values.  If 
these  two  steps  are  repeated  several  times,  an  arbitrarily  good  precision 
is  achieved  and  is  approximately  given  by  oy(Tm)//P,  where  P  is  the  number 
of  Ay  values  obtained  from  switching  back  and  forth.  Knowing  the 
character istics  of  both  the  random  and  the  systematic  deviations  of 
precision  oscillators  clearly  is  useful  to  the  designer,  the  planner,  the 
user,  as  well  as  the  vendor. 
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OSCILLATOR 


TEMPERATURE  SENSITIVITY  /  d*3  K 


Figure  2.  Nominal  values  for  the  tenperature  coefficient  for  the 
frequency  standards:  QU  -  quartz  crystal,  RB  =  rubidium  gas  cell, 
H  -  active  hydrogen  maser,  H(pas)  -  passive  hydrogen  maser,  and 
CS  -  cesium  beam. 


OSCILLATOR 


MAGNETIC  FIELD  SENSITIVITY  /  G 


Figure  3.  Nominal  values  for  the  magnetic  field  sensitivity  for  the 
frequency  standards:  QU  =  quartz  crystal,  RB  -  rubidium  gas  cell, 

H  “  active  hydrogen  maser,  H(Pas)  »  passive  hydrogen  maser,  and 
CS  =  cesium  beam. 
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OSCILLATOR 


REPRODUCIBILITY 


Figure  4.  Nominal  capability  of  a  frequency  standard  to  reproduce  the 
3ame  frequency  after  a  period  of  time  for  the  standards:  QU  =  quartz 
crystal,  RB  =  rubidium  gas  cell,  H  =  active  hydrogen  maser, 

H(pas)  ■=  passive  hydrogen  maser,  and  CS  =  cesium  beam. 


OSCILLATOR 


ABSOLUTE  ACCURACY 


Figure  5.  Nominal  capability  for  a  frequency  standard  to  produce  a 
frequency  determined  by  the  fundamental  constants  of  nature  for  the 
standards:  QU  »  quartz  crystal,  RB  =  rubidium  gas  cell,  H  *  active 
hydrogen  maser,  H(pas)  =  passive  hydrogen  maser,  and  CS  =  cesium  beam. 
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OSCILLATOR 


Figure  6.  Nominal  values  (ignoring  the  sign)  for  the  frequency  drift  for 
the  frequency  standards:  QU  =  quartz  crystal,  RB  =  rubidium  gas  cell, 

H  -  active  hydrogen  maser,  H(pas)  =  passive  hydrogen  maser,  and 
CS  *  cesium  beam. 
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Figure  8.  A  simulated  time  deviation  plot,  x(t),  with  indicated  ;sample 
time  t  over  which  each  adjacent  fractional  frequency,  y,,  is  measured. 

The  equations  are  for  the  standard  deviation  and  for  the  estimate  of  o  (t) 
for  a  finite  data  set  of  M  frequency  measurements.  It  is  often  the  case 
that  the  standard  deviation  diverges  as  the  data  length  increases  when 
measuring  the  long  term  frequency  stability  of  precision  oscillators, 
wehereas  oy(7)  converges. 
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Figure  11.  From  the  frequency  stability  characterization  shown 
in  Figure  10,  optimum  prediction  algorithms  to  minimize  the  time 
error  can  be  obtained.  Based  on  optimum  prediction  procedures 
the  RMS  time  prediction  error  for  a  prediction  interval  tp  can  be 
calculated  for  each  of  the  oscillators  shown  in  Figure  10,  and 
the  corresponding  values  are  plotted  in  Figure  11. 


Table  I. 


Typical  Applicable  Oscillators  and  Range  of  Applicabillt; 


a, 

Noise  Types 

Cs 

H-active 

H-passive 

Qu 

Rb 

2 

White  Noise  PM 

£  100  S 

Z  1ms 

1 

Flicker  Noise  PM 

$  1  s 

0 

White  Noise  FM 

l  10  3 

100  s  i-ti  104  s 

2  1  s 

21  3 

-1 

Flicker  Noise  FM 

£  days 

2104  s 

2  days 

2  1  s 

£10**  s 

-2 

Radora  Walk  FM 

2  weeks 

2  weeks 

i  weeks 

^  hours 

2  days 

Table  II 


Typical 

Noise 

Types 


Classical 
Standard 
Deviation 
of  x 


Classical 
Standard 
Deviation 
of  y 


White  Noise  PM 


Flicker  Noise  PM 


White  Noise  FM 


Flicker  Noise  FM 


t  •  o  (t)//3* 
(constant) 


t  •  oy(t) 


To  *  °y(To) 


Oy(T) 

~  °y(T) 


|2{N+1  ) 

V  3N 

2(N+1  Hnt 

i  3N 


°y( *0 ^ 


undefined 


,  /  N  In  N 
°yU  V  2(N-lHn2 

Oy(l)  /N/2 


Random  Walk  FM 


undefined 


Table  III 


Typical 

Noise 

Types 

Optimum 

Prediction 
x(x  ) 
rms 

Time  1 

Error  s 
Asymptotic 
form 

White  Noise  PM 

T  •  oy(-r)//3 

constant 

tint 

Flicker  Noise  PM 

-  ’  '  VT)  /list  1 

/£nTp 

White  Noise  FM 

t  *  Vt) 

V/2 

Flicker  Noise  FM 

t  •  Oy(t) /-/in? 

tp 

Random  Walk  FM 

T  *  Oy(t) 

-  3/2 
p 
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QUESTIONS  AND  ANSWERS 


PETER  KARTASCHOFF,  SWISS  PTT: 

Dave,  I  have  a  question  which  I  will  direct  to  you  and  to  Gernot  Winkler.  We  have 
measurements  and  we  have  recommended  designations  such  a  Y,  F  sub  Y,  sigma  sub  Y 
and  so  on.  Ten  years  ago  we  had  a  recommendation  in  the  CCIR,  International  Radio 
Consultative  Committee,  with  an  associated  report  which  explains  all  the  terms 
which  you  have  so  well  presented  here.  Everybody  uses  them.  Fifteen  years  ago 
there  was  a  subcommittee  of  the  IEEE  which  recommended  these  measurements.  Do 
you,  or  anyone  here,  know  if  there  is  any  other  published  international 
recommended  standard,  IEC  or  IEEE  besides  the  CCIR  recommendation?  As  far  as  I 
know,  the  CCIR  recommendation  still  stands  alone  and  I  do  not  know  what  the  IEEE 
is  doing  or  what  the  IEC  is  doing  about  this. 

MR.  ALLAN: 

As  far  as  I  know,  there  is  nothing  new  beyond  those  that  you  have  mentioned. 
There  is  an  IEEE  subcommittee  in  existence,  but  it  is  not  active.  I  know  of  no 
other  committee. 

A  COMMENT  HERE  BY  MR.  WINKLER,  NOT  INTO  THE  MICROPHONE  AND  NOT  DECIPHERABLE. 

MR.  BEARD: 

You  mentioned  environmental  effects  and  such  things  as  that  on  characterizing 
these  clocks.  Aren't  the  properties  of  the  measuring  system  with  which  you  make 
these  measurements  also  important? 

MR.  ALLAN: 

I  really  appreciate  your  mentioning  that.  Very  often  a  process,  in  fact  the 
millisecond  pulsar  is  an  example,  is  a  measurement  of  the  measurement  noise.  We 
think  that  we  cannot  get  the  real  millisecond  pulsar.  This  is  a  very  good  case  in 
point,  because  all  you  are  seeing  is  the  measurement  system.  In  many  applications 
that  can  be  a  problem.  This  goes  back  to  the  point  that  Dr.  Winkler  made:  one 
should  make  your  phase  or  time  measurement  system  as  good  as  possible  so  that  you 
don't  have  measurement  noise  contamination. 


